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Exact solution of the mixed-spin Ising model on a decorated square lattice with two 
different kinds of decorating spins on horizontal and vertical bonds 
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The mixed spin-(l/2, 5b, 5c) Ising model on a decorated square lattice with two different kinds 
of decorating spins 5b and 5c (5b 7^ 5c) placed on its horizontal and vertical bonds, respectively, is 
exactly solved by establishing a precise mapping relationship with the corresponding spin- 1/2 Ising 
model on an anisotropic square (rectangular) lattice. The effect of uniaxial single-ion anisotropy 
acting on both types of decorating spins 5b and 5c is examined in particular. If decorating spins 
5b and 5c are integer and half-odd-integer, respectively, or if the reverse is the case, the model 
under investigation displays a very peculiar critical behavior beared on the spontaneously ordered 
'quasi- ID' spin system, which appears as a result of the single-ion anisotropy strengthening. We have 
found convincing evidence that this remarkable spontaneous ordering virtually arises even though all 
integer-valued decorating spins tend towards their 'non-magnetic' spin state 5 = and the system 
becomes disordered only upon further increase of the single-ion anisotropy. The single-ion anisotropy 
parameter is also at an origin of various temperature dependences of the total magnetization when 
imposing the pure ferrimagnetic or the mixed ferro-ferrimagnetic character of the spin arrangement. 

PACS numbers: 05.50.+q, 05.70.Jk, 64.60.Cn, 75.10.Hk, 75.10.-b, 75.30.Kz, 75.40.Cx 
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I. INTRODUCTION 



One of the most fascinating and challenging topics 
in equilibrium statistical physics are phase transitions 
and critical phenomena of rigorously solvable interact- 
ing many-body systems [l[ . The planar Ising model rep- 
resents perhaps the simplest lattice-statistical model for 
which a complete exact solution is known since Onsager's 
famous solution 2] and the model simultaneously under- 
goes a non-trivial phase transition. Although the spin- 
1/2 Ising model on two-dimensional (2D) lattices with 
non-crossing bonds is in principle exactly soluble problem 
within the PfafRan method [3| , its precise treatment usu- 
ally requires an application of sophisticated mathemat- 
ical methods, which consequently lead to considerable 
difficulties when applying them to more complex models 
describing spin systems with interactions beyond near- 
est neighbors frustrated spin systems [f|, or higher- 
spin models with or without single-ion anisotropy and bi- 
quadratic interactions [6j] . It is noteworthy, however, that 
exact solutions of the spin- 1/2 Ising model have already 
been reported on several Archimedean lattices (square 
0, triangular and honeycomb 0], kagome extended 
kagome |9], bathroom-tile [Io| . orthogonal-dimer [ll[ and 
ruby [12J lattice) and more recently, also on a variety of 
more complex irregular lattices such as union jack [13j |. 
pentagonal [lij ]. square-kagome [lj|, as well as, two topo- 
logically different square-hexagonal [l6[ lattices. An im- 
portance of the rigorously solved 2D Ising lattices can be 
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viewed in providing guidance on the scaling and univer- 
sality hypotheses 17f and moreover, these exact results 
provide very valuable information about an accuracy of 
different approximative theories used to study spin sys- 
tems where the rigorous treatment is inapplicable. Even 
though the Ising model has been originally designed for 
describing essential features of insulating magnetic ma- 
terials [18], throughout the years various modifications 
of this model have found rich applications in seemingly 
diverse research areas (l9j . 

Over the last few decades, the mixed-spin Ising mod- 
els have attracted a great deal of research interest on 
behalf of much richer critical behavior they display com- 
pared with their single-spin counterparts. Actually, the 
mixed-spin Ising models are often convenient candidates 
for displaying tricritical phenomenon or other compli- 
cated change of an usual critical point to a multicritical 
point. Despite much efforts, there are only few exactly 
solvable mixed-spin Ising models, yet. Using the gen- 
eralized decoration-iteration and star-triangle mapping 
transformations, the mixed spin-(l/2, S) Ising model on 
the honeycomb, diced and several decorated planar lat- 
tices have exactly been treated by Fisher and Yamada 
many years ago [2fJ. Notice furthermore that an exten- 
sion of the generalized mapping transformations enabled 
to account also for the single-ion anisotropy effect. The 
effect of uniaxial and biaxial single-ion anisotropics has 
been subsequently exactly examined in the mixed-spin 
Ising models on three-coordinated honeycomb [2l| and 
bathroom-tile [22j lattices, diced lattice [23[ and some 
decorated lattices [2o| . To the best of our knowledge, 
these are the only mixed-spin planar Ising models with 
generally known exact solutions except several mixed- 
spin Ising models on the Bethe lattices studied within 
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the framework of exact recursion relations [26[. Among 
the striking models for which a precise solution is re- 
stricted to a certain subspace of interaction parameters 
only, one should further mention the mixed-spin Ising 
model on the union jack (centered square) lattice, which 
can be mapped onto a symmetric (zero-field) eight-vertex 
model with continuously varying critical exponents [27j ■ 

Exactly solvable Ising models on 2D decorated lattices, 
the bonds of which are decorated in various fashion by ad- 
ditional spins, are therefore of particular research interest 
(see Ref. [28| and references cited therein) . In the class of 
exactly solved decorated Ising models belong the origi- 
nal ferrimagnetic model introduced by Syozi and Nakano 
(29| , partly [30| and multiply 31 1 decorated models show- 
ing reentrant phase transitions, ANNNI models [32|], di- 
luted models of ferromagnetism [33[, decorated models 
with classical ^-dimensional vector spins [34J, Fisher's 
superexchange models and its numerous variants f35l , as 
well as, the models with higher decorating spins [201 ] . It 
is worthwhile to remark that rigorous solutions of these 
models have furnished answers to questions interesting 
both from the academic point of view (scaling and uni- 
versality hypotheses, reentrant phase transitions), as well 
as, from the experimental point of view (dilution, techno- 
logical applications of ferrimagnets). The vast potential 
of ferrimagnets with respect to technological applications 
has also stimulated exploration of the effect of single-ion 
anisotropy upon ferrimagnetic features of the mixed spin- 
(1/2, S) Ising models on the wholly [3] and partly [25| 
decorated lattices. 

Recently, Kaneyoshi [36| has proposed another notable 
example of the decorated Ising model on a square lat- 
tice the horizontal and vertical bonds of which are oc- 
cupied by two different kinds of decorating spins. Up 
to now, this remarkable model system have been stud- 
ied using the approach based on the differential opera- 
tor technique |37| . whereas an accuracy of obtained re- 
sults determined the Bethe-Peierls- Weiss approximation 
used for the undecorated lattice [36[. Therefore, the main 
aim of the present work is to extend the class of exactly 
solved Ising models by providing an accurate solution for 
this mixed spin-(l/2, Sb, Sc) Ising model on the square 
lattice with two different kinds of decorating spins Sb 
and Sc (Sb Sc). Exact results for the system under 
consideration are obtained by applying the generalized 
decoration-iteration transformation [28j establishing an 
exact mapping correspondence with an effective spin-1/2 
Ising model on the anisotropic square (rectangular) lat- 
tice. Owing to this fact, the known exact solution of the 
latter spin-1/2 Ising model on the rectangular lattice 
straightforwardly enables to acquire exact results for the 
former decorated mixed-spin model. Within the frame- 
work of this exact method, we will concentrate our atten- 
tion first of all to the influence of the single- ion anisotropy 
on the critical behavior and phase diagrams. Besides, 
temperature dependences of the total magnetization will 
be also particularly examined. 

The outline of this paper is as follows. In Section HT1 
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FIG. 1: A cross-section of the three sublattice (ternary) 
mixed-spin Ising model on the decorated square lattice. Open 
circles denote lattice positions of the spin-1/2 atoms (sublat- 
tice A), while black and hatched circles label lattice positions 
of the decorating spin-SB atoms (sublattice B) and the spin- 
Sc atoms (sublattice C), respectively. 



the detailed description of the considered model system 
is presented in the first instance. Then, some details 
of the mapping method will be clarified along with the 
derivation of exact expressions for the magnetization and 
critical temperatures. The most interesting results are 
presented and detailed discussed in Section [Hi] for two 
particular sets of the quantum spin numbers (5a, Sb, 
S c ) = (1/2, 1, 3/2) and (1/2, 2, 3/2). Finally, some 
concluding remarks are mentioned in Section llVl 



II. MODEL AND METHOD 

Suppose the square lattice with two different kinds of 
decorating spins Sb and Sc placed on its horizontal and 
vertical bonds, respectively, as it is diagrammatically de- 
picted in FIG.[T] As it can be clearly seen, the displayed 
magnetic structure can be also viewed as the three sub- 
lattice (ternary) mixed-spin system in which each vertex 
of the original square lattice is occupied by the spin-SA 
atom (sublattice A), while each its horizontal and verti- 
cal bond is decorated by the spin-SB (sublattice B) and 
spin-Sc (sublattice C) atom, respectively. As a result, 
the nearest-neighbor Ising model defined upon the un- 
derlying lattice is given by the Hamiltonian 
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2 ' 

c + 1, ■ ■ ■ , Sc denote three different kinds 
of Ising spins differing in value, their superscript de- 
termines the sublattice to which spins belong, whereas 
their subscript specifies lattice points where the spins 
are placed. The first two summations are carried out 
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over the nearest-neighbor A-B and A-C spin pairs, re- 
spectively, and the last two summations are taken over 
all decorating B and C sites. Accordingly, the param- 
eters Jab and Jac label pairwise exchange interactions 
between the nearest-neighbor A-B and A-C spin pairs, 
while the single-spin parameters Db and Dc measure a 
strength of the uniaxial single-ion anisotropy acting on 
decorating spins Sb and Sc, respectively. The essence 
of both single-ion anisotropy parameters Db and Dc lies 
in an uniaxial magnetic anisotropy, which comes from a 
crystal field of ligands from a local neighborhood of mag- 
netic centers [381 ] - 

The crucial step in our approach represents evalua- 
tion of the partition function, i.e. a sum over all possible 
spin configurations, what usually represents a formidable 
mathematical task in a highly cooperative spin system 
with many microscopical degrees of freedom. Fortu- 
nately, the problem connected with calculation of the 
partition function of decorated models can be simply 
avoided by the use of well-known trick, which consists in 
performing the summation over spin degrees of freedom 
of all decorating spins Sb and Sc before summing over 
spin degrees of freedom of the vertex spins S A . In doing 
so, the partition function can be partially factorized into 
the following product 



JV 



N 



z = e n^Bosft + s$) nwcish + s£ 2 ), (2) 



k=l 



where the symbol X){s A } stands for the summation over 
all available spin configurations on the sublattice A and 
the first (second) product is performed over all decorating 
spins (Sc) occupying the horizontal (vertical) bonds. 
Furthermore, the function Wq(x), which depends on two 
vertex spins S^ and Sfe from the sublattice A coupled 
indirectly via ith decorating spin, marks the expression 

+S S1 

Wq(x)= J2 exp({3D n n 2 )cosh(pj An nx), (3) 
n=— So 

where f2 = B or C, (3 = I/^bT), &b is being Boltzmann's 
constant and T is the absolute temperature. 

In order to proceed further with calculation, the gener- 
alized decoration-iteration mapping transformation [28| 
can be now employed 

Wa(S& + S&) = AnexptfRaS&Sb), fi = BorC. (4) 



transformation. As a matter of fact, the direct substitu- 
tion of four possible spin combinations of two vertex spins 
Sf^ and Stf indeed gives just two independent equations 
from the formula ([¥]), which subsequently unambiguously 
determine the mapping parameters An and Rq 



An = [Wn(l)W n (0)]* 
/3Rn = 21n[Wh(l)/Wh( 



= BorC, (5) 
ft = BorC. (6) 



At this stage, let us substitute the transformation (j4|) 
into the expression ^ in order to gain the relation 

Z((3,J A B,J A c,DB,Dc)^A^A^Z (f3RB,l3Rc), (7) 

which relates the partition function Z of the mixed- 
spin Ising model on the decorated lattice with the parti- 
tion function Zq of the undecorated spin-1/2 Ising model 
on the anisotropic square (rectangular) lattice with two 
different nearest-neighbor couplings Rb and Rc in the 
horizontal and vertical directions, respectively. Never- 
theless, it should be mentioned here that the mapping 
relation ([7]) between both the partition functions rep- 
resents a central result of our calculation, since it for- 
mally completes an exact solution of the partition func- 
tion Z with regard to the known exact result for the 
partition function Zn of the spin-1/2 Ising model on the 
rectangular lattice [2| given by the effective Hamiltonian 
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Now, we turn to calculation of the sublattice magne- 
tization. With the help of commonly used mapping the- 
orems [39|, one easily proves a validity of the following 
relation for the spontaneous magnetization m A of the 
sublattice A 



m A = (Si) = (Si ) = to , 



(8) 



where the symbols (...) and (. . .) denote canonical en- 
semble averaging performed within the decorated and its 
corresponding undecorated model system, respectively. 
Apparently, the sublattice magnetization m A directly 
equals to the magnetization mo of the corresponding 
spin-1/2 Ising model on the rectangular lattice with the 
effective horizontal and vertical coupling constants given 
by the equation ©. When taking into account the 
Potts's and Chang's [40j exact result for the spontaneous 
magnetization of the spin-1/2 Ising model on the rectan- 
gular lattice, it is possible to write the following expres- 
sion for the sublattice magnetization m A 



The physical meaning of the mapping Q is to remove 
all interaction parameters associated with one decorat- 
ing spin (5b or Sc) and to replace them by a new ef- 
fective interaction (Rb or Rc) between the remaining 
vertex spins fi-f and S^, . It is worthwhile to remark that TOA 
a self-consistency condition of the mapping relation (|4|) 
unambiguously determines both unknown mapping pa- 
rameters An and i?o, since it must hold independently 
of the spin states of both vertex spins included in this 
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or equivalently, 
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16W3(l)Wg(0)Wg(l)Wg(Q) 



[Wi(l) 



Wg(0)]* 



(9) 



(10) 



Next, it is also of particular interest to derive exact ex- 
pressions for the spontaneous magnetization of the sub- 
lattice B and C constituted by the higher decorating spins 
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5b and 5c, respectively. For this purpose, the exact 
Callen-Suzuki [41| spin identity formulated generally for 
both kinds of decorating spins can be successfully utilized 



(S?) = 



J2 S? Spexpj-PHj) 
J2s a exp(-/3H l ) 



(11) 



in which = B or C and the Hamiltonian Hi labels such 
a part of the total Hamiltonian |T]), which contains all 
interaction terms of one decorating spin 5p. By mak- 
ing use of the mapping relation ([7]) and the exact spin 
identity (fTTj) , it can be easily proved that the sublattice 
magnetization tob and toc satisfy the equality 



n = BorC, 



(12) 



where the function Fq(x) is defined as follows 



F Q (x) 



22 71 exp(/3L»nn 2 ) sinh(/3 Jaqtix) 

n=— So 



Sn 



(13) 



2J exp(/3£»nn 2 ) cosh(/3J A nnx) 

n= — Ssi 



By employing the differential operator technique |37| 
together with the exact van der Waerden identity 
exp(a5 A ) = cosh(a/2) + 25, A sinh(a/2) into the equation 
(TT2")) . one readily finds that both the sublattice magneti- 
zation tob and toc can be expressed solely as a function 
of the sublattice magnetization toa, in fact, 

m n = <exp[(5 A + Sfi)d/dx])F n (x) = 2to a F (1). (14) 

The above result completes an exact solution for the sub- 
lattice magnetization tob and toc in view of the formerly 
derived exact result © [or (jTUJ)] for the sublattice mag- 
netization TO-A- 

Finally, let us turn our attention to the critical behav- 
ior of the model under investigation. It is quite obvi- 
ous from the equation (| 14[) that all sublattice magnetiza- 
tion tend necessarily to zero if and only if the sublattice 
magnetization toa goes to zero. Accordingly, the critical 
temperature can be located from the condition which is 
consistent with the Onsager's critical condition for the 
spin-1/2 Ising model on the rectangular lattice Q 



sinh(/3 CJ R B /2) sinh(/3 CJ R c /2) = 1, 



(15) 



where /3 C = \/(k^T c ) is being defined as the inverse crit- 
ical temperature in energy units, T c denotes the critical 
temperature and the effective coupling constants Rb and 
Rc are given by the equation ([6]). It is also interesting 
to mention that the critical temperature can be obtained 
from the alternate condition, which is of course equiva- 
lent to the condition ([15]) . but without referring to the 
effective coupling constants of the undecorated rectangu- 
lar lattice. The expression (TIT))) for the sublattice magne- 
tization toa indeed yields the following critical condition 

{[^b(1)] 2 - [^b(0)] 2 }{[^c(1)1 2 ~ [^c(O)] 2 } 

= 4Wl(l)W§(0)WS(l)Wg(0). (16) 



The superscript in the aforementioned expressions means 
that the inverse critical temperature j3 c enters into the 
relevant expressions © instead of (3. 



III. RESULTS AND DISCUSSION 

Before proceeding to a discussion of the most inter- 
esting results, let us make few remarks on a validity of 
the results to be obtained in the preceding section. No- 
tice initially that the obtained results are rather general, 
actually, they hold regardless of whether ferromagnetic 
or antiferromagnetic interactions Jab and Jac are con- 
sidered, irrespective of the value of the decorating spins 
(5b and 5c) and even both the single-ion anisotropy pa- 
rameters Db and Dq can be taken independently of each 
other. It is noteworthy, however, that several studies 
reported on previously have already involved some par- 
ticular cases of the model under investigation. By im- 
posing 5b = 5c, Jab = Jac and Db = -De, for instance, 
our results reduce to those acquired for the mixed-spin 
Ising model on a symmetrically decorated square lattice 
adapted to study essential features of the ferrimagnetism 
[24]]. In the present paper, we will therefore restrict our 
attention only to the particular case with two different 
kinds of decorating spins 5b ^ 5c- More specifically, 
one of the two decorating spins (say 5b) is assumed to 
be integer valued, while the other one (say 5c) is an- 
ticipated to be half-odd-intcgcr. For simplicity, another 
constraint introduced through an equality between the 
single-ion anisotropy terms Db = Dq = D will be sup- 
posed in order to reduce the number of free parameters 
involved in the model Hamiltonian |T]). Other particu- 
lar models in which both the decorating spins are integer 
or half-odd-integer, respectivel y, w ill be explored in the 
subsequent separate works [42l. |43| . 

First, let us take a closer look at the ground-state 
behavior. In FIG. [H the ground-state phase diagrams 
in the D/Jab — Jac /Jab plane are depicted for two 
particular models with (5a, 5b, 5c) = (1/2,1,3/2) and 
(5a,5 b ,5c) = (1/2,2,3/2). As one can see, the inter- 
action parameters Jab and Jac energetically favor the 
highest spin state of the decorating spins 5b and 5c, 
respectively. Contrary to this, the easy-plane single- 
ion anisotropy {D < 0) has a tendency to lower the 
spin states of the decorating atoms. An eventual spin 
arrangement is thus determined by a mutual competi- 
tion between the exchange interactions and the single-ion 
anisotropy. It is quite obvious from FIG. [5] that there is 
considerable similarity between both the displayed phase 
diagrams. As a matter of fact, the only difference consists 
in the phase boundary that accompanies the transition 
2 «-> 1 of the spin-5B atoms, which is of course miss- 
ing in the first phase diagram. Finally, it is worthy to 
note the general condition allocating ground-state phase 

5 Q -1 



boundaries associated with the spin change 5 



D 



s si ,s n -i 



1 



1 - 25^ ' 



Q = BorC. 



(17) 



5 



f= (3/2,1/2) 
4 f 



• (order-disorder) 




(a) 



D/J. 



GO 



D/J. 



0.5 \ 0.0 

(-1/2,2,1/2) 



FIG. 2: Ground-state phase diagrams in the D/Jab — Jac/Jab plane for two particular spin cases with (Sa, 5b, Sc) = 
(1/2,1,3/2) [FIG. Ufa)] and (Sa,Sb,Sc) = (1/2,2,3/2) [FIG. Wib)]. Round brackets indicate spin order to emerge within 
different sectors of the phase diagram of the pure ferrimagnetic model with Jab > and Jac > 0. Different symbols to 
be defined in the legend characterize a change in the spin state that occurs at the displayed phase boundaries. Note that a 
sign change in Jab (Jac) leads merely to a sign change in the spin state of the spin-SB (Sc) atoms. The hollow-square line 
represents transitions between the ordered and disordered phases (for details see the text). 



Next, let us investigate in detail the effect of the single- 
ion anisotropy upon the critical behavior of the consid- 
ered model system. For this purpose, we have plotted in 
FIG. [3] the finite-temperature phase diagram in a form of 
the critical temperature vs. single-ion anisotropy depen- 
dence for the two investigated spin systems. It should 
be pointed out that the unique solution of the critical 
condition (fT5|) [or equivalently |[T6|) ] gives all the dis- 
played phase boundaries, which consequently represent 
the lines of second-order (continuous) transitions between 
the spontaneously ordered phases and disordered (para- 
magnetic) phase. Furthermore, it is worth noticing that 
the interaction parameters Jab and Jac enter within the 
critical conditions (fT5|) and (fT6|) only into the arguments 
of even functions. In this regard, the depicted phase 
boundaries remain in force irrespective of a sign of the 
interaction parameters Jab and Jac> what means that 
the same phase boundaries hold regardless of whether 
the pure ferrimagnetic, the mixed ferro- ferrimagnetic, or 
even the pure ferromagnetic system is considered. 

The most striking finding resulting from the finite- 
temperature phase diagram shown in FIG. [3] is that the 
lines of phase transitions do not reach zero tempera- 
ture at the boundary value of the single-ion anisotropy 
D\fil Jab = — 1, which is sufficiently strong to bring all 
the integer-valued decorating spins 5b into their 'non- 
magnetic' spin state S B — 0. Namely, one would intu- 
itively expect that all the phases with the 'non-magnetic' 
decorating spins Sb (i-e. the phases appearing in FIG. [2] 
in the parameter space where D < -D^o) should be in fact 
disordered due to their 'quasi-lD' character. In the zero- 
temperature limit, one actually finds that the effective 



coupling constant (3Rb substituting the integer-valued 
decorating spins Sb tends to zero, however, the effective 
coupling constant (3Rc substituting the half-odd-integer 
decorating spins Sc tends to infinity. In the consequence 
of that, the 2D decorated mixed-spin system indeed be- 
haves as the 'quasi-lD' spin system, since it effectively 
splits into a set of the independent mixed spin-(S , A, Sc) 
Ising chains. On the other hand, one should bear in mind 
that the system is spontaneously long-range ordered if 
and only if sinh(/3i? B /2) sinh(/3i?c/2) > 1, while it be- 
comes disordered just as sinh(/3i?B/2) sinh(/3i?c/2) < 1 
in accord with the critical condition (|15p . In the limit of 
zero temperature, it can be easily verified that 



lim 



sinh 



'0Rb 



I sinh(^y^ 



_ r oo if D > -Jab - \ Jac 
" I if D< -Jab - § Jac 



what means that the threshold single-ion anisotropy be- 
low which the system becomes disordered (at zero as 
well as any non-zero temperature) is given by the con- 
straint -D -d = — Jab — \ Jac- This non-trivial phase 
boundary, which cannot be obtained from simple ener- 
getic arguments, is depicted in the ground-state phase 
diagram (FIG. [2} as a hollow-square line. This result 
straightforwardly proves that all the phases appearing 
above this order-disorder line are at sufficiently low tem- 
peratures spontaneously long-range ordered in spite of 
their 'quasi-lD' nature. In agreement with the afore- 
mentioned arguments, the non-zero critical temperatures 
to be observed in the parameter space where D < D\ q 
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FIG. 3: Critical temperature as a function of the single-ion anisotropy for several values of the ratio |Jac|/|Jab| in two 
investigated spin systems with (S A ,S B ,Sc) = (1/2,1,3/2) [FIG. EJa)] and {S A ,S B ,S C ) = (1/2,2,3/2) [FIG. Mb)]. Different 
symbols characterize the same spin transition as explained in the legend of the ground-state phase diagram shown in FIG. [2] 
Broken lines connecting different spin transitions are guides for eyes only. 



manifest the order-disorder transition between the spon- 
taneously ordered phase and the disordered phase even 
though all the spin-SB atoms reside in the ground state of 
the ordered phase the 'non-magnetic' state S B = (see 
FIG. [3]). In this respect, the part of second-order phase 
transition lines starting at D = — Jab and terminating at 
D = - Jab - 5 Jac (cf. FIGs. [3 and ED can be identified 
as the critical line at which the spontaneous long-range 
order of the 'quasi-lD' spin systems disappears. 

In order to provide an independent check of the afore- 
mentioned scenario, the total and sublattice magnetiza- 
tion will be detailed analyzed for the decorated model 
with (Sa,S b ,S c ) = (1/2,1,3/2). Before starting our 
further analysis, it is worthwhile to mention that this spin 
system exhibits all generic features of the models with 
other decorating spin values, as well. The total magneti- 
zation mx = 1 77i a + 777b + mc| reduced per one atom of 
the sublattice A is shown in FIG. [4] for the pure ferrimag- 
netic system with Jab > 0, Jac > [FIG. Ufa)] and the 
mixed ferro-ferrimagnetic system with Jab > 0, Jac < 
[FIG. 1Kb)]. Note that different line styles are used to 
distinguish between different spin orderings to emerge in 
the zero-temperature limit (ground state). Since the sign 
change in Jac leads merely to the sign change in the sub- 
lattice magnetization mc (mc is an odd function of Jac)j 
all temperature variations of the total magnetization can 
be explained via the same temperature variations of the 
sublattice magnetization. 

For this reason, we depict in FIG. O all three sublat- 
tice magnetization as a function of the temperature for 
those particular values of the single-ion anisotropy for 
which the most significant thermal variations of the total 
magnetization occur. It is quite apparent from FIG.[5Ja) 
that a steep decrease in the total magnetization of the 



ferrimagnetic system to be observed at D/ Jab = —0.45 
can be explained through the vigorous thermal excita- 
tions | — > ~ of the decorating spins Sc, which can be 
clearly seen in the temperature dependence of the sublat- 
tice magnetization mc- On the other hand, the opposite 
thermal excitations \ — ► I are responsible for a rapid 
increase of the sublattice magnetization mc [see for in- 
stance the case D/Jab = —0.55 shown in FIG. EJb)], 
which in turn causes a gradual temperature-induced in- 
crease of the total magnetization. As far as the mixed 
ferro-ferrimagnetic system with Jab > and Jac < is 
concerned, the initial decrease in the total magnetization 
at D/Jab = —0.45 can be again attributed to the vig- 
orous thermal excitations | — ► | of the decorating spins 
Sc- However, it is easy to observe from FIG. G3a) that 
the thermal excitations of the decorating spins Sb over- 
whelm in the region of moderate temperatures (i.e. the 
sublattice magnetization tub declines more rapidly with 
increasing temperature than the sublattice magnetiza- 
tion mc), what consequently leads to a slight increase of 
the total magnetization in the mixed ferro-ferrimagnetic 
system (remember that the sublattice magnetization ttib 
is now oriented in opposite to the sublattice magnetiza- 
tion 771a and mc)- At D/Jab — —0.55, the total mag- 
netization of the ferro-ferrimagnetic system starts from 
zero because the sublattice magnetization me effectively 
cancels out both the sublattice magnetization m A and 
mc- The observed temperature- induced increase of the 
total magnetization can be mainly related to the ther- 
mal excitations | — > | of the decorating spins Sc, which 
are reflected in the temperature-induced increase of the 
sublattice magnetization mc- 

Another notable temperature dependences of the to- 
tal magnetization occur in the vicinity of the boundary 
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FIG. 4: Temperature dependences of the total magnetization for the decorated mixed-spin system with (Sa, Sb, Sc) = 
(1/2, 1, 3/2) when the ratio | Jac\/Jab = 1-0 is fixed and the single-ion anisotropy strength D/ Jab varies. FIG.|4ja) shows the 
situation in the pure ferrimagnetic system with the exchange constants Jab > and Jac > 0, whereas FIG. HJb) corresponds 
to the mixed ferro-ferrimagnetic system with Jab > and Jac < 0. Different line styles correspond to different phases. 




FIG. 5: The absolute value of the sublattice magnetization mA (dotted lines), uib (solid lines), and mc (dashed lines) for 
|Jac|/Jab = 1 and several values of the single-ion anisotropy parameter D/Jab = —0.45 [FIG.G^a)], —0.55 [FIG-Ob)], —0.9 
[FIG. 0c)], and —1.1 [FIG.[5fd)]. 
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value -Di,o/Jab = —1 at which the decorating spins 5b 
change their spin state. As a result, the relatively sharp 
decrease of the sublattice magnetization tob to be ob- 
served at D/Jab = —0.9 leads in the ferrimagnetic sys- 
tem to the relevant decrease of the total magnetization 
[FIG. Ufa)], while in the mixed ferro-ferrimagnetic sys- 
tem [FIG. BJb)] it is responsible for an effective increase 
of the total magnetization since the sublattice magneti- 
zation toa and mc are almost completely independent 
of temperature [FIG. 03c)]. Last but not least, the most 
interesting temperature dependences of the total mag- 
netization could be expected in the parameter space, 
where the 2D decorated mixed-spin system should effec- 
tively behave as the 'quasi-lD' system. FIG. EJd) shows 
the sublattice magnetization for the particular value of 
single-ion anisotropy D / Jab = — 1.1, which is sufficiently 
strong to force all the integer-valued decorating spins Sb 
towards their 'non-magnetic' state S B = at zero tem- 
perature. As one can see from this figure, the sublattice 
magnetization tob indeed starts from zero in agreement 
with our expectations and it is merely the effect of tem- 
perature that causes a steady rise of the sublattice mag- 
netization tob on behalf of promoting the — > 1 exci- 
tations of the decorating spins Sb- Altogether, it might 
be concluded that the shape of temperature dependence 
of the total magnetization is in the ferrimagnetic sys- 
tem almost entirely determined by the sublattice magne- 
tization tob (the sublattice magnetization toa and mc 
almost completely cancel out), while in the mixed ferro- 
ferrimagnetic system the rising magnetization tob lowers 
the mutual contribution of the sublattice magnetization 
toa and toc to the total magnetization. 



IV. CONCLUDING REMARKS 

In the present article, the generalized decoration- 
iteration mapping transformation has been utilized to 
obtain an exact solution of the mixed spin-(SA, Sb, 
Sq) Ising model on the decorated square lattice with 
the two different kinds of decorating spins 5b and Sc 
placed on its horizontal and vertical bonds, respectively. 
Within the framework of this exact method, the sought 
exact solution for the model under investigation has 
been attained by establishing a simple mapping relation- 
ship with the corresponding spin-1/2 Ising model on the 
anisotropic square (rectangular) lattice whose exact so- 
lution is known since Onsager's pioneering work (2j. In 
addition, the applied transformation method is of im- 
mense practical importance, since this method is rather 
general and it enables further interesting extensions. Ac- 
tually, it is quite straightforward to extend the applied 
procedure to account for: (i) the interaction between the 
next-nearest-neighboring spins Sa; (h) the multi-spin in- 
teraction between the decorating spins and their nearest 
neighbors; (iii) the biaxial single-ion anisotropy acting 
on the decorating sites; (iv) other decorated lattices such 
as decorated honeycomb or triangular lattices with (two 



or even three) different decorating spins; (v) decorated 
lattices with two or more decorating spins per one bond; 
etc. It is noteworthy, moreover, that the applied method 
is in principle applicable also to 3D decorated lattices, 
but unfortunately we cannot present an exact solution of 
the 3D decorated models because of the unknown exact 
solution of the corresponding spin-1/2 Ising model on the 
3D lattice (however, some results with a high numerical 
accuracy are available even for 3D lattices [44|). 

The most interesting result to emerge from the present 
study consists in providing an exact evidence for the 
spontaneous long-range order, which surprisingly appears 
in the 2D decorated spin system in spite of its 'quasi-lD' 
character. As a matter of fact, we have found a con- 
vincing evidence that the 2D decorated spin system re- 
mains spontaneously long-range ordered even if all the 
integer-valued decorating spins Sb are driven by a suffi- 
ciently strong (but not too strong) single-ion anisotropy 
towards their 'non-magnetic' state S B = and the sys- 
tem becomes 'quasi- ID' due to the effective splitting into 
a set of independent mixed spin-(S , A, Sq) chains. This 
finding has obvious relevance to the understanding of 
the 'quasi-lD' spin systems prone to spontaneous long- 
range ordering below some critical temperature, which 
necessarily need not arise from interactions establishing 
3D connectivity (3D magnetic lattice), but it may repre- 
sent an inherent feature of the 'quasi-lD' system. It is 
worthwhile to remark that this outstanding feature can- 
not be found in the decorated spin system with both the 
half-odd-integer decorating spins Sb and Sc [43], while 
it might be found in the decorated spin system with both 
the integer-valued decorating spins Sb and Sq provided 
that they are coupled to the vertex spins Sa through two 
different exchange interactions Jab ^ Jac, respectively 
[42| . From this point of view, the decorated spin system 
with the integer and half-odd-integer decorating spins Sb 
and Sq exhibits the spontaneous ordering related to an 
appearance of the 'quasi-lD' spin system the most nat- 
urally and hence, this exactly solved model would be of 
major importance if some experimental realization of it 
would confirm an appearance of the spontaneous ordering 
(and consequently non-trivial criticality) notwithstand- 
ing of its 'quasi-lD' character. 

Although the Ising model description may not be fully 
realistic for true magnetic materials, it is quite reason- 
able to expect that the exact solution of this simpli- 
fied model illustrates many important vestiges of the 
real critical behavior. Besides, the exact solution may 
also bring other valuable insights into the thermodynam- 
ical properties (magnetization, entropy, specific heat) of 
real magnetic materials without being affected by any 
crude and/or uncontrollable approximations. According 
to this, the main stimulus for the study of the mixed 
spin-(S , A, Sb, Sc) Ising model on the anisotropically dec- 
orated square lattice can be viewed in connection with 
possible experimental realization of this remarkable and 
rather curious magnetic lattice. It is therefore of par- 
ticular interest to mention that polymeric compounds, 
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which have the decorated square network assembly as the 
magnetic lattice, have been rather frequently prepared 
in an attempt to design novel bimetallic coordination 
compounds. The magnetic structure of the decorated 
square lattice is indeed peculiarity of two numerous series 
of polymeric coordination compounds with the general 
formula [Ni(L) 2 ] 2 [Fe(CN) 6 ]X.nH 2 [H and A[M B (L)] 2 
M A (CN) 6 ].nH 2 (M A = Fe, Mn, Cr, Co; M B = Mn, Fe) 
46| . In the former series, the magnetic Fe 3+ (Sa = 1/2) 
ions reside the square lattice sites and Ni 2+ (Sb = 1) ions 
decorate each its bond, while in the latter series the high- 
spin M B + ions like Mn 3 + (S B = 2) or Fe 3 + (S B = 5/2) oc- 
cupy the decorating sites and the low-spin M A + ions such 
as Fe 3+ (Sa = 1/2), Mn 3+ (S A = 1), or Cr 3 + (Sa = 3/2) 
reside the square lattice sites. Even though we are not 
aware of any trimctallic polymeric compound whose net- 
work assembly would consist of two different kinds of 
magnetic metal ions (decorating spins Sb ^ Sc) placed 



on the square net made up by the third magnetic metal 
ion (vertex spin Sa), the vast number of the bimetallic 
coordination compounds from the aforementioned series 
gives us hope that a targeted synthesis of such trimctallic 
compound (tailored from two structural derivatives with 
different decorating magnetic ions) could be successfully 
accomplished in the near future. 
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